Superconductivity Fluctuations in a One-dimensional Two Band Electron-Phonon 

Model with Strong Repulsive Interactions. 
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We study a one-dimensional, two band model with short range electron-electron repulsions (onsite 
U and nearest-neighbour V terms) and electron-phonon coupling. We show that there is a region of 
U,V and band filling in which singlet superconductivity fluctuations are dominant. This region is 
absent without electron-phonon interactions and includes large values of U,V. 
PACS numbers: 71.30 + h, 74.10 + v, 74.20 - z 



The physics of low dimensional strongly correlated 
fermion systems with repulsive interactions is a topic of 
active interest, largely due to the absence of a clear un- 
derstanding of the origin of high-T c superconductivity in 
cuprate oxides and the role of phonons in these correlated 
systems. In the framework of a simple one-dimensional 
(1-d) Cu-0 chain model @, we investigate the effects of 
both short range electron-electron (e-e) repulsive interac- 
tions (onsite U and nearest neighbour Cu-0 repulsion V) 
and electron-phonon (e-p) coupling on the ground state 
of the system. We show that superconducting (SC) cor- 
relations are absent in the model if we take into account 
e-e interactions only. The inclusion of e-p interactions 
leads to the appearance of a region of (U, V,p) (p, the 
band - filling) in which superconducting fluctuations are 
dominant. On the other hand the ground state of the 
system in the absence of e-e repulsion is a charge-density 
-wave (or spin-density-wave) state without a divergent 
SC response. Thus, the origin of the region with dom- 
inant SC response is the combined effect of c-e and e-p 
interactions for this model. 

We use a renormalization group (RG) two cut-off ap- 
proach, developed in earlier works With some as- 
sumptions on the model parameters our analysis is valid 
in the large U, V limit. The possibility of superconduct- 
ing fluctuations in quasi- 1-d systems with strong repul- 
sive e-e interactions and e-p coupling was first raised in 
work of Zimanyi et.al. where results are obtained for 
a massive Thirring model. The two-band model without 
e-p interactions was considered in jjj, where numerical 
results are presented, pointing out the existence of SC 
fluctuations in a strong coupling limit. 

We consider a chain consisting of two types of atoms: 
Cu on odd sites with d-orbitals and O on even lattice 
sites with p-orbitals. The Hamiltonian of the system is 



H = Hq + H e 



(1) 
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where t is the hopping integral, < i,j > are nearest- 
neighbour sites, A = (E p ~ Ed)/2, E p and Ed are site 
energies, Ud, U p are Hubbard onsite repulsive energies, 
and V is the repulsion amplitude between nearest neigh- 
bour sites. Direct antifcrromagnetic coupling betveen Cu 
sites is omitted. Also 



Hep — H e p \ -\- H e p 2, 



(4) 



and we consider two models of electron-phonon coupling: 
the molecular crystal (MC) model with the Hamiltonian 
H ep i in which optical phonons couple to the electron site 
energy; and the Su-Schrieffer-Heeger (SSH) model with 
the Hamiltonian H ep ,2 in which the lattice distortions 
modulate the electron-hopping matrix element t. The 
Hamiltonian H ePt i consists of two parts H ep i = H ep d + 
H e p,p, where the each part has the form 



H ep — 
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Kq~ + Xqipi = 
(d k +d^)p k , 
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with Wo = (k/M) 1 / 2 , g = \/(2Mu;o) 1/2 , Pk = E4+ q c q - 
(M ion mass, loq optic phonon frequency, k elasticity 
constant, A e-p coupling constant.) All terms in @ have 
indices d or p and the sum is over odd or even sites for 
Hep.di H e p^p respectively. The Hamiltonian H ep ^ takes 
into account intermolecular phonon modes 

P 2 1 

H ePt2 = E ?if + o K (*+i - ~ E Sti 



<1,3> 



(2) 



where Sti t j = X(qi — qj), uj q — 2(k/M) 1 / 2 sin(ga/2) is 
the acoustic phonon frequency, 



(G) 



1 



g(k,q) = i4Asin(ga/2) cos(fca + qa/2)/ (2Mui q ) 1 / 2 , and 
a is the Cu-0 lattice constant. 

First we consider the noninteracting Hamiltonian Ho- 
Diagonalization gives 



H = V^t 2 cos 2 ka + A 2 (c+(/c)c 2 (fc) - Cl (fc)+ci(fc)), 



where 



Cd(k) = cos^feCi(fc) 4- sm9kC 2 {k) 
c p (k) = — sin6>fcCi(fc) + cos9kC 2 (k) 



(7) 



(8) 



with tan(20 fe ) = -2t cos(fco)/A, -7r/2 < 28 k < tt/2. We 
now have a 2-band electronic structure and consider the 
case of an entirely filled lower band. The filling factor 
of the upper band is < p < 2 ( empty for p = and 
filled for p = 2). With unit cell 2a, the quasimomen- 
tums k and k + ix /a are equivalent, and we may con- 
sider that states in the lower band have quasimomcnta 
in the interval —n/2a < k < -n /2a and in the upper band 
n/2a < \k\ < 7r/a; then k F a = tt/2 + 7rp/4. The Fermi 
velocity is 



v F 



2at sm(2k F a) 
^J\t 2 cos 2 (k F a) + A 2 ' 



(9) 



Since we will use an RG approach we take into account 
only states in the upper band in the vicinity of Ep which 
are described by operators c%. Then Ho has the form, in 
x-representation, 

H Q = V F y+ + (-i — )* 2 , + + VpV+_(i — )* 2 ,-, (10) 

where ^ 2i ± includes momenta near ik F , respectively. 
Below we will omit the index 2 and also terms in H 
comprising (taking into account terms with \&i can 
produce a shift of the chemical potential and some renor- 
malization of the Fermi velocity.) Therefore we can sub- 
stitute in the Hamiltonian terms 

^ d — > sm9 F ^(x), * p — > cos9 F (11) 

Note that in the cases t/A < 1 or 1 we have 

isin(7rp/4) 



sin#; 



(12) 



We first consider e-e interaction effects. For the Cu- 
O case it is appropriate to consider Ud ~3> U p . We first 
study the case U p — 0. The effect of small U p is easily 
taken into account and will be discussed below. In terms 
of a "g-ology" model J3J the Hamiltonian H ee gives the 
scattering amplitudes 



9i 



Ua 



sin 4 9 F + 2Va sin 2 9p cos(2kpa) = 53 



92 = — — sin 4 9p + 2Va sin 2 9p 



94, 



(13) 



where gi is the back scattering amplitude, and 32,34 are 
forward scattering. The "Umklapp" part 33 exists only 
for the half filled case p = 1: for the sake of simplicity we 
will not consider this case. Since we use a RG approach 
below, we consider gi/ixvp < 1, that is Ua,Va < irvp 
or sin 9 <C 1 for large U, V. We have spin - rotation 
invariance, i.e., g± = g^. Therefore, if it is not essential 
we are omitting spin indices. The effect of the 174 term is 
taken into account separately: it simply produces a shift 
in the velocity of the spin and charge degrees of freedom: 
v a =v F (l + gi), v p = v F (l - g^ 

The usual RG equations defining the scaling behaviour 
of the system are H 



9'i 



1 



9\ 



g c — gi — 2g2 = const. 



(14) 



(15) 



For g\ > the excitation spectrum is gapless g\ — > 
g\ = 0, and there is a gap if g\ < 0. The charge ex- 
citation spectrum is gapless if g c > and has a gap 
A p if g c < 0. The ground state has a most divergent 
singlet (triplet) superconductivity (SC) response when 
9c > 0, gi < (<7 C > 0, g\ > 0). In our case 

g c = sin 4 9 F + 2Vasm 2 9 F cos 2 9 F (cos(2k F ) - 2) < 0. 

(16) 

Therefore there is no region (U, V) with divergent SC 
fluctuations. Possible ground states are charge- or spin- 
density wave, depending on the sign g\. (This sign can 
vary due to the cosk F a term). We see that in order 
to obtain SC correlations it is necessary to have large 
positive gl or negative g% terms. As we will see below 
this condition can be achieved by taking into account e-p 
interactions. 

Second order perturbation theory in e-p interaction 
produces a retarded e-e interaction j^j for u> less than 
a Debye frequency w < uj d ~ (k/M) 1 / 2 . (We con- 
sider the case ujjj < Ep). The effective e-e interac- 
tion can be described in "g-ology" terminology: gi, p h = 
-2g 2 (k F ,2k F )/u 2 k F , g2, P h = -2g 2 (kp,0)/uo, gz, P h = 
gi.ph (half- filled band only). In the case of the MC model 



(17) 



(18) 



The parameters k, A in p7|) , (p.8|) are, of course, differ- 
ent, as well as other similar parameters in Hamiltonians 



gi,ph — 9%ph — g3,ph — 

whereas the SSH model (^) gives 

A 2 

9\,ph = 93, P h = -4 — (sin 2 9 F cos 2 9 F ). 

K 
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H ep i, H ePt 2- Note that all terms are negative, and g^'ph 
is due solely to onsite e-p coupling and does not contain 
renormalization terms sin#i?, cos9p. In the case 9 < 1 
the onsite e-p interaction is dominant. We now have two 
types of e-e interactions with cutoffs Ep and cup,- Wc 
thus use the RG procedure for a two cut-off model. 
The one-loop scaling equations (|l4|) , ([l5|) for the <?j are 
unaffected by the presence of retarded interactions. The 
equations for the gi, p h, taking into account cross terms 
9i9j, P h, were derived in ||: 



9i, P h 



1 



f 



.3 

-(ofli + o3c + gi,ph)gi, P h 

7TV p Z Z 



52, ph 



0. 



(19) 



(20) 



We will consider the case g^^h = 0. The integration 
in (p^[),(|20|) is taken from Ep to uj q ~ wd(wo), where 
^d(^o) is the renormalized value of Lop, Q. As a result 
the combined action of different scattering processes is 
described by 



9i 



9i 



9i : ph ' 



(21) 



The properties of the system at energies small compared 
to luq are derived from a model with single interactions 
gj and bandwith ojq. 

We now examine the solutions of equations (|l4|), (15), 
( |l9| ) , (po| ). The initial conditions for (|l^), (|l5|) are defined 
by (|13|). The initial conditions for (|l9|),(|2C|) are defined 
by We denote g { ^ h = -7, g® h = -7. If 

.9i°' > (we will see that this is the situation in the 
interesting region), from (|l^) we obtain that g\ scales 
toward small positive values g\ <C gf^ . Note that in the 
case 9 < 1 we have gf ph ~ g^ph^ that 1S 7 ~ 7- F rom 
(|l|) it follows that g{ - 2g* = g { ° ] ~ 2 5 ° 2) . A positive 
derivative in (|l9|) implies that gi tP h scales toward large 
negative values. We consider the opposite case g[ ph < 0. 
Then, at least initially, gi tP h will scale toward a small 
negative value. Therefore we demand that 



~ 2 9 X 



(°) j. i (o) 

2 9c 



A>o, 



(22) 



since g\ lP h < 0. The inequality (|22|) can not be valid 
throughout the scaling process, since g\ scales to small 
values. Therefore the value g\ ph may not be very small. 
We do not require \g\ ph \ <§; 7; for our purposes it is suf- 
ficient that gl ph > 7 — 27, as we show below. The value 
92,ph is not scaled as follows from (p0|), i.e., g^ ph = — 7. 
This value does not contain the renormalization coeffi- 
cient sin^F- As a result of scaling we have the state with 
gf = g* + g* ph - The ground state of the system with 
the new scaling amplitudes has dominant divergent SC 
susceptibility if 



9c 



T 

9i 



2ff 2 T 



91 



(0) 



2<?f + 9*i,p h + 27 > 0. (23) 



Since we suppose that gf ph « g\ h < 0, we have a state 
with spin gap A a . Therefore the dominant singularity is 
the singlet SC response with SC correlation function 



R(x) 



' l+gj/2nv p 
1 - gT/2nv p 



> 1. 



(24) 



In this case the CDW response can be divergent with cor- 
relation function oc x~ Kp . The inequalities (|22|) , (f23|) de- 
fine the region in which singlet SC correlations are dom- 
inant. In terms of u = ^ sin 4 9, v = Va sin 2 9 cos 2 9 we 
rewrite (||),(||) as 

7 + 2t;(l + 2cos(-y)) < u < 2 7* - 2u(2 + cos( 

(25) 

where 27* = 2j + g^ ph . It is easy to obtain the solution 
of (H). This is the region ABCD in Fig.l delineated by 



lines: u ~ 0, v = 0, « = 7 — 2v, u = 27 + g\ 



ph 



2v. 




C D 



FIG. 1. Phase diagram obtained by two cut-off RG scaling. 
Divergent SC response regions are I and II. 

Recall that in the case 9 <C 1, Eq. (|l2|) , the bare repul- 
sive energies U ~ u/9 4 , V ~ v/9 2 ^ 7. Thus our model 
includes the case of strong electron repulsion. For any 
point (u, v) in the region ABCD (p5[) is valid for 



p > — cos 

7T 



(max{- 



7 — 2v 27* — 4u — u 



4u 



2v 



}, (26) 
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In the case f/A <C 1 we can obtain the phase dia- 
gram in term of the bare values U, V. Then the coor- 
dinates of points A,B,C,D are A = {0, (47*)(A/t) 4 }, 
B = {0,2 7 (A/t) 4 }, C = {( 7 /2)(A/t) 2 ,0}, D = 
{( 7 */2)(A/t) 2 , 0}. The SC region is deformed to include 
region II due to the sin(7T/o/4) term. The equation of the 
curve EF is 

v A 2 ((2 7 *- 7 )fc + 4 7 + 16 7 *) 2 9 , T/ A 2 
V= ^ 72(fc + 2)( 7 + 2 7 *) ' = T 2 "' (2?) 

In the limit k —* oo we have f oc V 2 , but in this region 
p ~ 1/J7 1 / 4 -> 0. The inequality (§|) with t/A < 1 
becomes 



(U 1 + 8V l )y 2 -6V 1 y-j>0, 
{AV X - U)y 2 - 6V lV + 2 7 * > 0, 

where Vi = V(A/t) 2 ,[/i = J£(A/t) 4 , y = 
In region I the solution of (ES) is po < p 
sin 2 (7rpo/4) = yo is the largest root of eqn. 
region II we have p\(U, V)< P < P2{U, V), where pi, pi 
are easy obtained from ( |28|) . If V = the solution is 



(28) 

sin 2 (7rp/4). 
< 2, where 
(M). In the 



1/4 



1/4 



for {Ji > 7; if J7 = , then in the region 7 /2 < V\ < 7 * 
the solution is sin 2 (-7r/9/4) > yo, where yo is the largest 
root of eqn. @ for U = 0. 

In using the RG approach, we supposed as usual 
that gi/nvF < !• For small t/A we have the ini- 



(0) 



(t/A) sin(7rp/2). Recalling that <7i 



tial value 

(t sin(7rp/4)/A) 2 or (t sin(7rp/4)/A) 4 , gi, v h = const, we 
suppose that our results are reasonable if we are not too 
close to band edges, that isei<p<2 — £2, and p ^ 1 
(.93 = 0). It follows from our analysis that in region III 
we have large spin and charge gaps, so that there is only 
a charge-density-wave divergent response. In the region 
IV we have g^ < , small gf < and thus divergent 
charge- and spin-density-wave (if g\ — -» 0) responses. 

In our consideration above we have not taken into ac- 
count effects of Up repulsion. This is easily achieved by 
substituting u into ( |2p| ) in the form u = ([/a sin 4 6p + 
U p a cos 4 9f)/2- For small values of U p , the RG approach 
remains valid, and all results in terms of the new u, v are 
retained. For t/A 1 we have cos#f ~ 1, so that we 
can not consider the large U p limit in our approach. 

In conclusion, using a two cutoff RG approach we have 
studied a two band, 1-d tight-binding model with e-e 
and e-p interactions. We included onsite U and near- 
est neighbour V electron repulsions, as well as intra- and 
inter-molecular e-p coupling. We have shown that there 
is no U, V, p, t, A parameter region with dominant diver- 
gent SC response in the absence of e-p interactions. In 
the lowest order RG approach we found that such a re- 
gion does occur if we include e-p coupling with optical 



intra-molecular modes. Only this form of e-p interaction 
produces an effective renormalized g2,ph term. We have 
found that the singlet SC region includes large values 
of the U, V repulsive interactions if tsin(7r / o/4)/A <C 1. 
Note that this behavior is impossible in the framework 
of a one-band model, for which A = 0. Then, instead of 
we have 



7 + 2^(1 - 2 cos(vrp)) < U < 2 7 + g* l ph - 2^(2 - cos(irp) 

(30) 

where 7 = -g [ ^ ph , 7 = -g^, < p < 2. The solu- 
tion of ( |30|) is the same region ABCD in Fig.l, provided 
2 7 + g\p h > 7 . However the bare values U, V must be 
very small, of the order of phonon scattering strengths. 
Moreover, in this case the requirement 2 7 + gl h > 7 
can be not realized. (Recall that in the two band case we 
have 7 ~ 7 in the limit t/d <C 1, since effective intersite 
phonon coupling terms are small due to the factor (t/A) 2 
). Note also that we have used a RG approach. There- 
fore we did not consider the strong coupling limit where 
a phase separation instability could take place our 
results are valid in some vicinity of the SC fluctuation 
phase. 

We have proposed one possible scenario for the ori- 
gin of dominant SC fluctuations in quasi-one-dimensional 
systems as a result of the combined effect of repulsive e-e 
and attractive e-p interactions in a 2-band situation. We 
suggest that features of this picture will survive in anal- 
ogous two-dimensional models of high-T c superconduc- 
tors [||, in particular in 3-band Peierls-Hubbard models. 
However the orbital structure of the order parameter in 
this case (s-wave vs d-wave) is unclear without detailed 
calculations. 
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